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Abstract 

We reconstruct a rational Lax matrix of size R+1 from its spectral curve (the desingularization 
of the characteristic polynomial) and some additional data. Using a twisted Cauchy-like kernel (a 
bi-differential of bi-weight (1 — v, u)) we provide a residue-formula for the entries of the Lax matrix in 
terms of bases of dual differentials of weights u, l—u respectively. All objects are described in the most 
explicit terms using Theta functions. Via a sequence of "elementary twists", we construct sequences 
of Lax matrices sharing the same spectral curve and polar structure and related by conjugations by 
rational matrices. 

Particular choices of elementary twists lead to construction of sequences of Lax matrices related 
to finite-band recurrence relations (i.e. difference operators) sharing the same shape. Recurrences 
of this kind are satisfied by several types of orthogonal and biorthogonal polynomials. The relevance 
of formulae obtained to the study of the large degree asymptotics for these polynomials is indicated. 



^Work supported in part by the Natural Sciences and Engineering Research Council of Canada (NSERC), Grant. No. 
261229-03 and by the Fonds FCAR du Quebec No. 88353. 
■^bertola@crm.umontreal.ca 

^Work supported in part by the NSF grant # 0400484. 
* mgekhtmaOnd . edu 



Contents 



1 Introduction and setting 1 

2 Notation and main tools 5 

2.1 Theta functions 5 

2.2 Prime form 6 

3 Inverse spectral problem for rational Lax matrices 7 

3.1 Using different tensor weights 9 

3.1.1 Twisting by flat line bundles 11 

3.2 Residue formulas for the Lax matrix 12 

3.2.1 Left eigenvector 13 

3.2.2 Structure of the Lax matrix near a branch-point of X 13 

3.2.3 Change of divisor 14 

3.2.4 Change of line bundle 15 

3.2.5 Linear (smooth) deformations of the line bundle Xri 15 

3.3 Spectral bidifferential 16 

3.4 Elementary twisting lattice and dual wave functions 17 

3.4.1 Dualization 19 

3.5 Expressions in terms of Theta functions 20 

4 Finite band recurrence relations and commuting (pseudo) difference operators 22 

4.1 Lax and ladder matrices 25 

4.2 Riemann-Hilbert problems 27 



5 Commuting (pseudo)-difference operators in duality related to the two-matrix model 28 
1 Introduction and setting 

The aim of this paper is to present an explicit solution of the inverse spectral problem for Lax matrices 
A{x) of size {R + 1) X [R + 1) depending rationally on x. The forward problem and beautiful 
connections with integrable systems were explored in [3, 14, 27, 28, 1, 2, 25, 21, 22, 23, 24, 16] (to 
name a few); in particular it was shown in these works the important role of the theory of Theta 
functions [15] in the solution of both forward and inverse problems. 

From the literature cited above we know that the forward problem (under suitable genericity 
assumptions) produces "spectral data" consisting of a smooth algebraic curve C of genus g with 
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two meromorphic functions X,Y : C ^ C, where X has degree R + I; X{p) = x is the spectral 
parameter while y = Y{p) is the eigenvalue of the rational matrix to be reconstructed. 

The scheme of reconstruction requires that we fix two dual tensor weights i^,! — u; this 
means that the eigenvectors of A{x) will be realized as bases in the suitable spaces of sections of 
z^-differentials (for the left eigenvectors) and (1 — i/)-differentials for the right eigenvectors. The 
parameter u can be chosen integer or half-integer. This type of Baker-Akhiezer functions was 
considered (but in a different context) in [29] and slightly earlier in the series of papers [22, 23, 24]. 
In addition we need a divisor T (whose degree depends on u) and an arbitrary meromorphic differential 
r]; they will determine the local properties of the dual BA vectors by fixing the zeroes/poles and the 
essential singularity structure, the latter determined by exp2i7r J rj. 

The problem is not new and Baker-Akhiezer functions have been around for a long time (see, e.g. 
surveys [13, 20] or the monograph [3]); however our aim is to provide explicit residue formulae 
for the entries of A{x) and explore the intimate relation between a suitable bidifferential ^(p, g) of 
weight (1 — u, v) and several objects of the theory. 

We then "twist" the reconstruction scheme so as to obtain recurrence relations for the BA vectors. 
We construct a sequence r„ of divisors of the same degree such that T^+i — r„ is an elementary 
divisor of degree zero consisting of two arbitrarily chosen points (in general). These elementary 
twists set in a general framework the original idea behind the construction of discrete variable BA 
functions suggested in the papers [21] and later utilized in [25] to develop the theory of commuting 
difference operators. 

The associated sequence of bidifferentials then defines a pair of sequences pn, i^n of v/{l — v)- 
differentials (called dual wave functions) related by a particular form of Serre duality: such 
duality is realized via a residue pairing of the form res PnT^m = ^mn, where oo^"*") is the divisor 

oo(+> 

of positive points in the elementary twisting divisors. 

If the elementary twisting divisors are chosen amongst the poles of X we obtain wave-functions 
solving a finite-term recurrence relation of the form 



These recurrence relations fall within the scope of the theory of "difference operators" extensively 
studied [27, 25]. Our interest has a different origin in connection with the theory of orthogonal 
and biorthogonal polynomials and their asymptotics for large degrees; in this perspective the wave- 
functions represent a (formal) asymptotic regime for the polynomials, the sequence of bidifferentials 
^ is intimately related to "Christoffel-Darboux-like" kernels and "Christoffel-Darboux-formulae" 
that arise in those contexts. 




(1-1) 
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We repeat the argument of [10] to illustrate this connection in the simplest case of the ordinary 
orthogonal polynomials. 

Orthogonal polynomials Pn{x) with respect to a weight on the real line, wn{x) = e~^^(^)dx, 
satisfy a three-term recurrence relation^ [30] 

XPn{x) = -fnPn+lix) + (3nPn{x) + 7„_ip„_i(x) , (1-2) 

which can be written in matrix form as xp = Qp, with p the semiinfinite vector of the orthogonal 
polynomials and Q the tri-diagonal (symmetric) matrix with entries given by the coefficients of the 
above recurrence relations. 

In studying the large degree asymptotics one typically sends the large parameter N appearing 
in the measure to infinity at the same rate as the degree n of the polynomial [12]: this means 
that -while Q implicitly changes because of the change in the measure- we are considering the 
polynomials and the recurrence relations very "far down" along the diagonal. On a heuristic level 
one argues that the tridiagonal semiinfinite matrix Q can be replaced by a doubly-infinite matrix 
X (i.e. indexed by Z rather than N) of the same shape and symmetries. Consider now the associated 
functions 7r„ := p„e~"2"^(^): if V{x) (the potent/a/ appearing in the measure that defines the OPs) is 
a polynomial of degree d+ 1 then this sequence -while still satisfying the same three-term recurrence 
relation- satisfies also a 2d + 1-term differential recurrence relation 

J^dxTTn = Cd{n)TTn+d + ■■■ + -Cd{n - (i)vr„_d (1-3) 

where in matrix form the recurrence is represented by a (skew-symmetric) matrix P with d supra- 
and d sub-diagonals. The scaling is needed (on heuristic grounds) to assure the boundedness of 
the coefficients of the recurrence relation. By construction, the two matrices P,Q satisfy 

\PM = ]^-^ (1-4) 

and in the N ^ oo limit they commute: we thus replace them by two commuting doubly-infinite 
matrices X, Y of the same shape and symmetries. At this point, the first problem is therefore to 
classify such pairs of commuting matrices and much of this has been extensively analyzed in [25]; 
some additional ingredients (Serre duality) can be found in [10] and are put in a general context in 
the present manuscript. 

In applications stemming from random matrices, the so-called Christoffel-Darboux kernel has 
crucial importance since it generates all correlation functions [26]. The C-D kernel is nothing but the 
orthogonal projection operator (for the chosen measure) on the subspace of polynomials of degree 
— 1 or less and is given by 

N-l 

Kn{x,x') = ^Pj{x)pj{x') . (1-5) 
j=0 

^In fact the situation allows a generalization to holomorphic weights on contours as explained in [10]. 
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Due to the Christoffel-Darboux theorem it can be expressed in terms of only two OP 

KN[x,x)=-fN (1-6) 

and this fact is crucial in proving universality results since it allows to express the asymptotic behavior 
for large N in terms of a fixed (i.e. N independent) number of polynomials (in this case 2). 

In the heuristic approach used in [10] (and then justified rigorously using Riemann-Hilbert tech- 
niques) the quasipolynomials 7r„ were replaced by meromorphic sections of a spinor bundle, namely 
by half-differentials on the (asymptotic) spectral curve, in this case hyperelliptic. The function x was 
then regarded as a meromorphic function X{p) on this algebraic curve whose multiplication of the 
half-differentials 7r„ can be expressed in term of the same sequence of half-differentials, thus produc- 
ing a recurrence relation. The "orthogonality" was replaced by a residue pairing between the sequence 
7r„ of half differentials and the Serre-dual sequence vr* of half-differentials : res7r„7r* = 6mn- Simi- 
larly the kernel Kn{x,x') was replaced by bidifferential of weights (1/2, 1/2) that played the role of 
projection operator with respect to the residue pairing. 

Such bidifferential also satisfies a "Christoffel-Darboux" theorem 

^ . vr7v(p)7r^„i(y) -^7V-i(y)^7v(p) 

^^(P^P ) = xip) - xip') ^'-^^ 

which is ultimately an identity for Theta functions; this is fully generalized presently in Prop. 3.2 
and Prop. 4.2. 

Our paper does not focus primarily on difference operators, rather we find them as a byproduct 
of the sequence of transformations induced on the Lax matrix by the elementary twisting; also, the 
eigenvectors for the Lax matrix (i.e. the Baker-Akhiezer vectors) solve certain Riemann-Hilbert prob- 
lems with quasi-permutation monodromies. These were studied in [19] for their own sake, while our 
approach finds them as a natural byproduct of the inverse-spectral reconstruction. Riemann-Hilbert 
problems with quasipermutation monodromies are also related to asymptotics of (multi)orthogonal 
polynomials; indeed after the so-called normalization of the RH problem satisfied by the polynomials 
and associated functions, one is lead to an approximating asymptotic problem with quasipermutation 
monodromies. 

The paper is organized as follows: in Section 2 we recall the basic tools from the geometry of 
Riemann surfaces, in particular the notion of Theta functions and prime forms, after [15]. 

In Section 3 we set up the inverse spectral problem for rational Lax matrices; here the problem 
is solved using pairs of dual Baker-Akhiezer vectors with tensor weights u,! — u where u € ^Z,. 
A residue formula for the Lax matrix in terms of spectral projectors is derived. We also derive the 
"generalized Toda lattice" in terms of elementary twists and express the ladder matrices and the 
matrices implementing the change of a line bundle in terms of suitable residue formulas. Finally, we 
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provide explicit expressions for the relevant twisted Cauchy kernels in terms of Theta functions and 
prime forms. 

In Section 4 we specialize the generalized Toda lattice so as to obtain genuine finite-terms re- 
currence relations (difference operators); in this setting more explicit formulas for the BA vectors 
are derived. The connection to Riemann-Hilbert problems with quasi-permutation monodromies is 
pointed out. 

Finally, in Section 5 we consider the case that is potentially most relevant to the study of biorthog- 
onal polynomials for the two-matrix model [6] and reveal a notion of duality that is well known for 
biorthogonal polynomials but was not known in the context of pairs of commuting difference opera- 
tors. 

We end this introduction pointing out that in the case that £ has genus all the formulas can 
be expressed in terms of rational functions of the uniformizing parameter: this is left as exercise for 
the interested reader. 



2 Notation and main tools 
2.1 Theta functions 

For a given smooth genus-gi curve L with a fixed choice of symplectic homology basis of a and 
6-cycles, we denote by the normalized basis of holomorphic differentials 



= ^jt , f ^1 = Tjg = Tij . (2-1) 

We will denote by Q the theta function 

e(z) := gi7rn-Tn-2i7rz-n ^2-2) 

For brevity we will often omit any symbolic reference to the Abel map: namely if p G £ is a point 
and it appears as argument of a Theta-function, it will be understood that the Abel map (with a 
certain basepoint) was applied. 

We denote by /C the vector of Riemann constants (also depending on the choice of the basepoint) 



9 r 

e=i 



rp J.: 



(2-3) 



where in this expression the cycles aj are realized as loops with basepoint po and the inner integration 
is done along a path lying in the canonical dissection of the surface along the chosen representatives 
of the basis in the homology of the curve. 
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The crucial property of /C is that for a nonspecial divisor T of degree g, T = X^j=i7j, the 
"function" 



/(p) = e(p-r-/c) 



(2-4) 



has zeroes precisely and only at p = ^j, j = I . . . g. 

We will also have to use Theta functions with (complex) characteristics; for any two complex 
vectors e,6 the theta function with these (half) characteristics is defined via 



e 



(z) := exp ( 2iTT { + . z + • 5 jj 6 I z + ^ + 



Here the (half) characteristics of a point C € O are defined by 

2C = (5 + re 



(2-5) 



(2-6) 



where the factor of 2 is purely conventional so that half integer characteristics have integer (half)- 
characteristics. This modified Theta function has the following periodicity property : for A,/^ € 



e 



z + X + Tfi) = exp iTT{e- \ — 6 ■ fi) — ivr/i ■ t ■ n — 2iTTz ■ fj, 



e 



(^) 



Note also the symmetry 



e 



(z) = e 



(-Z) 



(2-7) 
(2-8) 



2.2 Prime form 

The prime form E{(X') is defined as follows [15] 

Definition 2.1 The prime form E{C,C) is the (-1/2, -1/2) bi-difFerential on C x £. 

eA(u(C)-u(c')) 



/iA(C)/i[»](CO 



k=l 



(2-9) 



(2-10) 



u=0 



where ujk are the normalized Abelian holomorphic differentials, u is the corresponding Abel map and 



A 



is a half-integer odd characteristic (the prime form does not depend on which one). 



The prime form E{(,(') is antisymmetric in its arguments and it is a section of an appropriate line 
bundle, i.e. it is multiplicatively multivalued on C x C: 



E{C + a„ O = EiC, CO , EiC + b„C') = EiC, CO exp 



(2-11) 
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In our notation for the (half)-characteristics, the vectors a, (3 appearing in the definition of the prime 
form are actually integer valued. We also note for future reference that the half order differential 
is also multivalued according to 



3 Inverse spectral problem for rational Lax matrices 

The goal of this section is to explore the inverse spectral problem, namely how to reconstruct a 
matrix rationally dependent on x from the knowledge of its spectral curve and some additional data. 
Since the construction is quite symmetric we can also deal with the dual situation without any extra 
effort, thus treating the spectral parameter and the eigenvalues on the same footing. 
We work with the following data 

1. A smooth curve C of genus g. 

2. Two meromorphic functions X, Y with polar divisors X, 2) of degrees R + I and + 1. 

3. A (generic) divisors F of degree g + R (not necessarily positive). 

The main tool is the following adaptation of the Cauchy kernel [15, 22, 23, 24] 

Proposition 3.1 For a generic choice of divisor T there exists a unique kernel K{p,^) which is a 
function w.r.t. the point p and a differential w.r.t. the point ^ with the divisor properties 



such that res K{p, ^) = 1. The subscripts above indicate in which variable the divisor properties are 
considered. 

The proof follows easily from the Riemann-Roch theorem; we will write explicitly the expression of 
this kernel (in a generalized setting of which the current one is a particular case) in terms of Theta 
functions later on. 

Remark 3.1 [Linear equivalence] In fact we could be slightly more general in the formulation of the above proposition, 
since what matters there is only the equivalence class (modulo principal divisors) ofV — X. In particular we could use 
in (3-1, 3-2) two different divisors T>, —V (both of degree —g — 1) as long as V and T> are equivalent. 

In that case, however, res K{p,(,) = /(p) would be a meromorphic function with divisor (/) = T> — V; there is 



only one such function (generically) up to scalar multiplication. Hence the normalization would have to be fixed in 
some other ad hoc way. 



hA{p + aj) =e'^"^hA{p) 



(2-12) 
(2-13) 



{K{p,i)\>-T + X-i 
{K{p,Ok>r-X-p 



(3-1) 
(3-2) 
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Example 3.1 If the divisor X. consist ofR+1 distinct pointsxo, 
for K is 





PpOOQ iX) 




• Ppoooilg+R) 


Cdet 




Ppooh(7i) ■ 


■ Pp^nX^g+R.) 


^i(C) 


LJl(7l) 


^l(79 + ii) 




- ^AO 







. xe and V is a positive divisor, then the expression 



(3-3) 



where PpqiS) stands for the (unique) normalized Abelian differential of the third kind with first-order poles at p,q and 
residues ±1; the constant C depends on the divisors X, F and is chosen so that the residue at ^ = p is 1. 

Consider now the expression M{p,^) := {X{p) — X{^))K{p,^); its divisor properties w.r.t. p,^ 
follow from the properties of K: 

{M{p,0)p>-r (3-4) 

(M(p,e)5 >r-2X, (3-5) 

where the pole on the diagonal is now absent because of the multiplication by X{p) — X{^). 
Again the Riemann-Roch theorem implies that generically 

r{-T) = i{T)-g + degT + l = R + l (3-6) 

since (generically) i(r) = 0. 

Also i(r — 2X) = R + 1 (generically); to see this we note that the space of third-kind differentials 
with poles not exceeding 2X has dimension deg(2X) — I + g = 2R + 1 + g. Imposing the vanishing 
at g + R points gives as many linear constraints, hence reducing the dimension to i? + 1. 

Let i/jqIp), . . . ,tljfi{p) be any basis of the vector space of meromorphic functions with divisor 
exceeding —T and let <^o(p)i • • • > ^Fi.{p) be any basis of the vector space of differentials with divisor 
exceeding F — 2X. Let us introduce the notations 



'Poip) 
'^r{p) 



(3-7) 



These vectors will be called the pair of dual Baker— Akhiezer vectors and we will show later on 
how this term is motivated by the Serre duality. Note that the notion of the dual Baker- Akhiezer 
function was first introduced in [11] where it was applied to construct algebro-geometric solutions 
to Gelfand-Dickii, NLS and sine-Gordon hierarchies. 

Since the dual pair spans their respective spaces it follows that 



M{p, G C{^o, ■■■^r}® C{v9o, ...,^r} 



(3-8) 
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In other words there is a {R+ 1) x (R + 1) matrix IK with constant coefficients such that 

iX{p) - X{i))K{p, = 'P\0^'4^ip) (3-9) 

which gives immediately 

Proposition 3.2 There is a constant matrix IK of size R + 1, depending on the choice of bases 
ipj, ifj, such that 

From this expression we derive the following identity. 
Corollary 3.1 Independently of the choices of the bases we have the identity 

dX{p) = ^\p)^^{p) . (3-11) 
Moreover, ifp, q e C\{TUX) are two points such that X{p) = X{q) then (f\p)Kip{q) = 0. 
Proof. Taking the residue on the diagonal we have 

1 .esi.(M) = «M, ,3.12, 

5=p dX{p) 

The second statement follows from the fact that K{p,q) is regular and hence the numerator in its 
expression must vanish whenever X{p) = X{q). Q.E.D. 

The Lax matrix can now be constructed from these spectral data if tp and (p^K are viewed as the 
right/left eigenvectors with eigenvalue Y{p) at the point(s) p above x = X{p). The explicit residue 
formula will be given later in Sect. 3.2 in a generalized setting. 

3.1 Using different tensor weights 

In the above scheme we are using a pair of BA vectors with tensor weights and 1 respectively, that 
is, functions and differentials. This is, in fact unnecessary and in some applications (typically to the 
asymptotics of ODEs) it may even be too restrictive. 

In general, we could widen the scope of the construction so that i/> and (p can be tensors of 
weight u and 1 — u respectively; similar considerations (motivated by applications to quantum field 
theory) were used in [22, 23, 24, 29]. The tensor weight v can be typically integer or half-integer; the 
particularly useful case [10] is = ^ where both elements of the pair are spinors (half-differentials). 
Clearly some modifications in the way the Riemann-Roch theorem is applied are needed (i. p., more 
general Serre-duality arguments). 
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If hy{'D) is the dimension of the space of i^-differentials with divisor exceeding V then we know 
that, for v>\, K{V) = if P is generic and degP > 5yi + {2v - l){g - 1). 

Thus, let r be generic and of degree deg(r) = {2v — \){g — 1) + R + 1; then one finds from the 
Riemann-Roch theorem that, for v G ^Z, j > \ (the case v = \ being the one discussed above) 

/ii_^(-r) = K{T) + deg(r) -{2v-l){g-l)=R+l . (3-13) 

Thus we should use a divisor F of degree {2v — l){g — 1) + i? + 1 and choose t/j to be a basis in 
Hi^y{-T) (of dimension R + 1). 

The dual BA vector tp would then span H^iV — 2X) also of the same dimension. 

The relevant "Cauchy" kernel is then a bidifferential of weights (1 — v,v) with divisor properties 



{^{p,q))i.\ > ^ " ' (3-14) 




'{:} 

and normalized by the requirement that the 'V-residue" is one, namely the expansion along the 
diagonal p = q\n local coordinate z{p) = z, z{q) = z' is 

^{P, q) = , (1 + 0{z -z')). (3-15) 

z — z' 

The leading coefficient of the expansion is invariant under changes of a local coordinate and hence 
it is a geometrical quantity. 

Note that if is a half-integer then, for completeness, one should also choose a spinor bundle 
(i.e. signs with which the half-integer spinor changes along each handle). 

Completely similar considerations as before show that 

Proposition 3.3 There exists a constant matrix ¥S'^'> such that 
and hence 

dX{p) = ip{p)K^''^ip{p) (3-17) 

In the next sections we can now consider this more general case but remove the explicit reference 
to the tensor weight u. Thus we will have 

• Fa generic divisor of degree {2v — l){g — 1) + R + I; 

• = [tpQ^ . . . , tpji] a vector of basis (1 — i/)-d if ferentia Is in Hi_;^(— F) (of dimension R + 1); 

• (f = [ifQ, . . . , ifji] a vector of basis i/-differentials in H,y{T — 2X); 

• ^{p, q) the unique bi-tensor of bi-weight (1 — u, v) with the divisor properties and the normal- 
ization listed above. 
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3.1.1 Twisting by flat line bundles 

In applications to isospectral dynamics and in several other applications it is necessary to consider 
a slight generalization of the above picture, in that one twists the line bundle implicitly associated 
to the divisor F by some other line bundle (which may depend on external parameters or "times"). 
Typically (as in the classical example of finite-gap integration of KP or KdV dynamics) the extra 
line-bundle data falls within the class to be described below. 

Let T] be an arbitrary meromorphic differential on the curve C with pole divisor 

K 

{r]) >'^diCi, , ci,- ■ ■ ,CK e C , di, . . . ,dK eN. (3-18) 

i=l 

Let us denote its residues by U 

K 

:= resr/ , Vtj = . (3-19) 

Ci ^ — ' 

i=l 

The Abelian integral 2171 rj (where the base-point of integration affects only an overall normaliza- 
tion) has in general nontrivial periods around the 2g handles of the curve and around the punctures 
Ci. The exponential e^'^'^^" '^ defines a homomorphism of 7ri(£ \ {ci, . . . , c/^-}) i— > and hence a 
certain flat line-bundle. Moreover this line-bundle has transition functions of exponential type^ at 
the punctures q. 

Twisting the previous description by this line-bundle £^ is then equivalent to considering v- 

differentials (resp. (1 — i/)-differentials) with essential singularities at the punctures Cj, i = 1, . . . , K 

of the same type as e^'^'^'^ . 

The twisted i/-Cauchy kernel is then a (1 — z^, i/)-bidifferential with singularities of the form 
^2tnJ''ri g^(j g2«7r/^»?^ ^^^|.^ ^^^^ ^^jH same symbol) 

exp (^2i7r -q^ K{p, q) (3-20) 

is locally a (z^, 1 — z^) bidifferential with divisor > — F + X — q{\np) and >r — X — p{\nq), and with 
multiplicative multivaluedness along the homotopy group of £ \ {ci, . . . ck} given by the character 

:vri(£\{ci,...,cx}) ^ C 

X^(7) =exp(^2i7r jfr/^ . (3-21) 

and such that near Cj in a local coordinate z it has a singularity of type 2^^*'^*^. Thus in general, 
unless the residues ti are integers, this kernel has logarithmic branching at the points q. 



®If di = 1 then the singularity may be a pole or powcr-likc singularity with non-integer exponent and hence also branching 
singularity. 
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Correspondingly, the bases (ip) are 1 — i^-differentials (resp. (z>)-difFerentials) with essential 
singularities of type e^^*'^-^ The uniqueness of such kernel is a simple argument in function 
theory and Riemann-Roch theorem. On the existence we do not insist at this point (although it 
would not be difficult to prove it abstractly) since we are going to produce explicit expressions in 
terms of Theta functions in Sect. 3.5. 



3.2 Residue formulae for the Lax matrix 

Let Z be the polar divisor of Y. We start with the observation that Y{^)'ip{^)^{p, ^) is a 1-differential 
(in ^) with poles only at 5,X and a simple pole at p with residue —Y{p)rp{p), therefore 

Y{pMp) = -resy(e)^(e)J^(p,0 = Y] resy(e)V'(0-^(P,e) = 



iP{p) (3-22) 



The expression 

is -a priori- a rational expression in x\ it has poles at the X-projection of the divisor 2 of poles of 
Y and at x = cxo (i.e. has a polynomial part). 

In particular if t) is a pole of Y of order /c with oo ^ Xq = X{t}) and does not coincide with any 
branch-point of X (i.e. dX|j5^_i^^ ^ ^ 0) then A{x) has a pole of order /c 

A{x) = 0{{x - Xo)-^) (3-24) 

If t) is a branch-point of X and /^t > 1 is minimum order of branching of X (/^ = 2 being the case of 
a simple branch-points) then 

A{x) =o[{x- Xo)-[*^/'^l) (3-25) 

If t) coincides with one of the poles of X then the Lax matrix will have polynomial parts of degree 
k or [k/{d — 1)] (if d is the order of the pole). Obviously the degree of A depends on the maximal 
degree amongst all poles of Y above x = oo. \\ Y has no poles coinciding with any of the poles of 
X then A{x) will be necessarily bounded at x = oo. 
By construction we have 

Y{p)ii^{p) = A{X{p))-iP{p) (3-26) 

Therefore (as expected) ip{p) is the right eigenvector of the matrix A{X{p)) with eigenvalue Y{p); 
the different points p lying above the same values of X{p) give the (generically distinct) eigen- 
value/eigenvector pairs. 
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3.2.1 Left eigenvector 

Consider now Y{p)R{p,^)ip{py: this is a 1-difFerential in p with poles in pat^,^ and ^. Since the 
z^-residue at ^ = p of ^{p,^) is 1, it follows from a simple computation in a local coordinate that 
the residue at p = ^ of this bidifferential is Y{^)ip'^{^). Therefore 

Y{0^\0 = res r(p)^*(p)i?(p,e) = -Yl resy(p) y/?^t^f,V *(p) = 

= ^'(0 res = (3-27) 

AW:= E-i-W^^^™ (3-28) 

^-^ p=q X — X[p) 

It is clear from the defining formulaethat 

I(x)IfC = -KA{x) (3-29) 
Therefore the left eigenvector of A{x) is (^*(^)]fC : 

Y{Oip\0^ = ^\OMX{0)^ = ^\O^Ai^) ■ (3-30) 

From Corollary 3.1 it follows that 



is the normalized left eigenvector { a v — l-differential) of A{x), in the sense that 

l{p) ■ tPip) = 1. (3-32) 

From the second part of Corollary 3.1 follows also that (as it should) the evaluation of i{q) at the 
other points q ^ C above X{p) are orthogonal to i/'(p). Note that the dual left-eigenvector has poles 
at the branch-points of the X projection. 

The (generically) rank-one projector on the eigenspace with eigenvalue Y{p) is given by 

n(rt:=^(rt8<(p) = ^<5!0>5 (3-33) 
3.2.2 Structure of the Lax matrix near a branch-point of X 

Suppose c is a critical value of X and X^/^j^i = X~^{c). Let us choose local coordinates near the 
point as 

Zi = {X-c)^^ . (3-34) 
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Then we have , . , , , ^ . x 

We see that each residue extracts the fij jet of Y and 11, contributing to a rank-//j Jordan block of 
the Lax matrix and this point. 
Indeed, by expanding 

y (e) = yj,i^/m) = E n.-^^/d^. jes ^^^§^ = ^.,M.-i-^n,/ (3-36) 

If IV{z)dz = R{z)L{z)dz, one sees easily by induction that the rank of any linear combination of the 
first k derivatives of 11 is less or equal to k, and generically is precisely k. 

3.2.3 Change of divisor 

The matrix A{x) depends implicitly on the divisor F but its characteristic polynomial does not (since 
the latter describes the algebraic relation between the rational functions X,Y on the spectral curve 
£). We investigate how a change in the divisor T (within the same class of generic divisors) affects 
the Lax matrix ^4(2;). 

Proposition 3.4 Let Ar{x) the Lax matrix constructed as in Section 3.2 using a divisor T. Let F 
be another divisor of the same degree and with similar genericity properties; then 

A^{x) = C^^^{x)A^{x)C~l{x) (3-37) 

where 



C7rr(^)=£es '^:^^V.;r (3-38) 



^P^r(-)=^^^^^^ = ^r.r(-)-^ (3-39) 

Proof. The product 'i/j{£,)^{p,C) has poles in ^ only at T, X and at ^ = p with residue —ip{p). 
Hence -as before- _ 

V'(p) = res ^{miv, = res ^rV^^^^(p) (3-40) 

from which the expression for Cp f{x) follows. The expression for Cf p(x) and the fact that it is the 
inverse of Cpp(2;) follows from simply interchanging the roles of F and F. Analogous expressions 
can be found for the change of divisor for the matrix Ar{x) (using the basis of forms (p). Q.E.D. 

Note that the transition matrices C^f{x) are rational functions of x with poles only at the X- 
projection of the divisor F and with a polynomial part of degree equal to the degree of the subdivisor 
of F that coincides with some poles of X. 
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3.2.4 Change of line bundle 

The matrix A{x) also depends on the chosen third-kind differential t] and we investigate what 
happens when we vary the differential within the same class. 

Namely, let rj be another third kind differential and let ip, -0, K, K be the same objects constructed 
before but using the differential rj instead of t]. 

The one-form (in q) 

K{p,q)^{q) (3-41) 

has poles at p = q with residue —ip{p), at X and essential singularities at the poles of i],r]. The 
transition matrix is then given by deformation of contours 

^{p) = res ]K(p, q)^{q) = res tfY^'^^M p) = M,,^{X{p))^l,{p) (3-42) 
with _ 

Here by the symbol res we simply mean the integral around a small loop encircling the point (the 
differential is not meromorphic there but has essential singularities). The matrix jy(x) has thus 
essential singularities (in general) in the complex x-plane at the points X{cj),X(cj). 
Consequently the Lax matrices are related by a simple conjugation. 

A^{x) = M^^^{x)-^A^{x)Mr,,j^{x) (3-44) 

Note that -by construction- the Lax matrices are still rational with the same pole structure even 
if related by a conjugation with a non-rational matrix. These formulae provide the integration of 
any isospectral dynamics on rational matrix-valued functions available in the literature on integrable 
systems (see, e.g. [14, 1, 2, 28] and references therein) and are, in fact, even more general. Indeed 
if r] depends (smoothly) on one or several a "time" parameters, the above formulae would provide 
integration of the flow on the isospectral manifold induced by the dependence of 7]. 

Note that the framework we are proposing is more general than the one in [1, 2] since the 
differential t] need not have poles coinciding with any of the poles of X (which is the case in loc. 
cit.). 

3.2.5 Linear (smooth) deformations of the line bundle Xv 

If 7] depends linearly (or even smoothly) on a set of times generically denoted by t then Mj^M~^ is 
a rational matrix (i.e. without essential singularities) as long as the residues of rjt are independent 
of t. 
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To show this suppose r] = rjt depends smoothly on a parameter t; in the literature the dependence 
is taken to be linear in the sense that the coefficients of the singular parts at the poles (in some 
chosen and fixed local coordinate) evolve linearly in t but the statement we are making here is more 
general in that it may include any deformation, including a motion of the position of the poles. 

Let iptjift be the dual bases of sections evolving in a smooth way; the reader should realize that 
this evolution implies a "gauge arbitrariness" consisting in the freedom of (smooth) change of basis 
within the same vector spaces of sections. This arbitrariness makes no difference on the rational 
nature of the infinitesimal deformation. 

Denoting by a dot the "time" derivative we note that ijj (and (p) are then sections of the 
same tensor space with additional singularities. Indeed, near any of the poles cj of rj we have 
tp = ftiv)^^"^* , with ft{p) analytic near Cj 

^l^=( rvtft + f^e^"^' . (3-45) 



Thus 



Note that f^rjt near a pole of r]t has a pole singularity (without logarithmic term) because the 
residues of 7]t are independent of t and hence rjt is a second-kind Abelian differential. 
Applying the argument that have already used several times before, we obtain 

V'i(p) = - res ^t{p, q)tl^tiq) . (3-46) 

q=p 

The expression we are taking the residue of, is a differential with poles only along X + p and along 
the divisor of poles of rjt -which we denote by due to r]t. Thus 

*,(P) = - res A(p. ,)*,(,) = ,|a^ t'M-'x(> ''''>- (^"^ 

M,(.)M.(.r'= r« MM?)l£ (3-48) 

qeX+€ X - X{q) 

Since the differential in the numerator is meromorphic (without essential singularities) on the 
spectral curve £, the latter expression is rational in x E C. 

3.3 Spectral bidifferential 

Suppose we are given a rational {R + 1) x (R + 1) matrix A{x) and let us denote by ya{x) its 
(generically simple) R+1 eigenvalues; consider the following bidifferential 

S[[x,ya[x)y, [x ,yb{x ))) := ^ (3-49) 

{x-xY ^((^ _ yJ(3;))Tr((A - y,){x')) 
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where the tilde denotes the matrix of co-factors (the classical adjoint). Since 



n,(:c) := ^^ZJtl^"^^ (3-50) 
T,{{A-ya){x)) 

is the (generically rank-one) projection onto the eigenspace with eigenvalue ya{x), it is not difficult 
to see [9] that this bidifferential extends naturally to the spectral curve and, in fact, we are now going 
to prove this in general. This object has appeared in the context of isomonodromic deformations 
of rational connections in the sense of [17]; indeed in [9] it was shown that this bidifferential is the 
generating function of the Hessian of the logarithm of the isomonodromic tau function. 

Here we are not considering such deformations but we can relate easily this bidifferential to the 
i^-Cauchy kernel introduced above. Let p,q be the abstract points on C with coordinates {x,ya{x)) 
and {x' , yh{x')). Note that the expression (3-49) seems to have at first sight a double pole whenever 
two points project to the same x-value, but in fact this occurs only when the branch of the eigenvalue 
is the same, meaning that it is a double pole only on the diagonal of the symmetric product of the 
spectral curve with itself. 

Indeed it follows immediately that since Ila{x)dx = '4>{p)ip^{p)K (and letting x = X{p),x' = 

n) - Tr ( '^(p)v>'(p)^ ii^{q)^\q)K \ _ y*(p)IKV>(g) y*(g)lCV>(p) _ 
^^^'^^ - [x{p)-Xiq)X{p)-Xiq)) - Xip) - X{q) X{p) - X{q) " W^'^') 

(3-51) 

so that the spectral bidifferential is nothing but the symmetric square of the z^-Cauchy kernel. It will 
be shown that S{p,q) is the square of the Szego kernel in Cor. 3.3; notice that, for the time being, 
the symmetric square has only a double pole on the diagonal p = q and no other singularities (this 
follows from the divisor properties of K{p,q), (3-14) and the type of essential singularities. 

3.4 Elementary twisting lattice and dual wave functions 

Suppose that we specify a sequence of elementary divisors 'Xn of degree 0, n € Z; by "elementary" 
we mean that they are of the form 

T„ = cx>i+) - ooi") (3-52) 

where oon^\ ooi ^ are two sequences of points arbitrarily (but generically) chosen. We will assume 
that, for any m and re, oo^^ ^ ooi~^'s are distinct from each other (but points within the same 
sequence may be repeated). 

If we twist the "initial" divisor T 

Tn :=r 
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r„:= + (3-53) 

r„+i - T„+i n < -1 

we obtain a sequence of divisors r.„ of degree {2v — l){g — 1) -|- ii -I- 1. Same strategy as before can 
still be applied (generically), namely we have a corresponding sequence of bases ■0^ and (f^ and of 
Christoffel-Darboux kernels ^n{p,0 3" satisfying 

Let us define cxj^"*") as the set of all points ooi^^ (counted without multiplicity with which the 
may appear in our sequence). 

Definition 3.1 The divisor oo(+) = |J{ool^^} will be called the dualization divisor. 
Proposition 3.5 The sequence of kernels Rn{p,q) satisfies 

{0 TTi Tl 

^m{p, q) - ^n{p, q) m> n 

The proof is a simple inspection of the residues; the product of the two kernels is a differential in 
^ that has no poles on oo(+' if n > m; and if n < m then it has only poles at a finite number of 
points of oo(+) and at = p,^ = q, where it has the indicated residues. 

Consider now the difference Kn+i{p,q) — ^n{p,q)', since both kernels have z^-residue 1 on the 
diagonal, this difference is regular there. Moreover 

{^n+l{P, q) - ^n{P, q))p > - 00^^\ + X (3-56) 

(j^+1 {p, e) - J^„(p, q))q > r„ - ooi-\ - X (3-57) 

The two divisors on the right hand side of these inequalities have degree —{2v — l){g — 1) — 1 
and {2u — l){g — 1) — 1 respectively; it follows that (generically) there is a unique meromorphic 
(1 — i/)-differential 7r„ and a unique meromorphic i^-differential pn in the respective spaces specified 
by these divisors. We have proved that 

Corollary 3.2 The difference of two consecutive Christoffel-Darboux kernels in the generalized Toda 
sequence factors 

•^+i(p, q) - ^{p, q) = T^n{p)Pn{q), (3-58) 

with 'Knip): Pn{q) defined (up to multiplicative constants) by the requirements 

Mp)) > -r„ - ooi+\ + X (3-59) 
(/On(<?)) >r„-ooi~\-X (3-60) 



18 



By induction, 



n+L-l 
j=n 



(3-61) 



Remark 3.2 The name of "generalized Toda lattice" is due to the fact that if X is the projection of a hyperelliptic 
curve (with two simple poles C 3 oo± above x = oo) and we use the sequence of elementary divisors = oo+ — oo_ 
then we recover a setting of the standard Toda lattice theory by looking at suitable isospectral evolution. 

Using the Christoffel-Darboux kernels one can therefore reconstruct a sequence of Lax matrices 
An{x) all sharing the same spectral curve and connected by conjugation by the transition matrices 
Cn,m{x) := C'r„,r„(a^) introduced in Prop. 3.4; these ladder matrices satisfy the obvious relations 
(which entail discrete integrability) 

Cn,m{x)Cm,e{x) = Cni{x) , Vn^m^ieZ. (3-62) 

Cn,UX{p))lPn{p) = ^miP) (3-63) 
ipi{p)KnCnAX{p)) = 'pUP^m (3-64) 



Remark 3.3 The previous construction is a generalization of the "discrete variable Baker-Akhiezer function", an 
idea originally formulated in [21] which is the hinge of the theory of commuting difference operators. 

3.4.1 Dualization 
Note that 

res PmTTn = Smn (3-65) 



oo' 



(+) 



since if m 7^ n the product is a differential with a polar divisor of degree 2 supported only at the 
points ooli^^'s or only at the points oon ^'s; only if m = n one pole is at oor^'^ and one at oo^ ' so 
that the residue over all 00^+^ is nonzero. The fact that this residue is actually 1 follows from Prop. 
3.5 

res (j^+i(p,0 - ^n(p,0) (■*^m+i(C,g) - ^{^,q)) = Smn {iin+i{p,q) - ^{p,q)) (3-66) 

which implies that res PnTr^ = 1. 
In addition we have 

res pm{p)iin{p,q) = I ^ - ^.^^ 7rm{q)Sin{p,q) = I ~ (3-67) 

oo<+) pm m < n 00'+) U m < n 
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3.5 Expressions in terms of Theta functions 



We now present explicit expressions of all the objects introduced so far. Let us decompose the divisor 
r (of degree (2i/ - l){g - 1) + R + 1) as 



2u-l 



r = 71 + . . . + + Yl ^0 



(3-68) 



where Fq are divisors of degree g — 1 and 71, . . . ,7_r+i are R points singled out arbitrarily. 

Let ^1, . . . C2u~i be arbitrary fixed points (the final formulae will only have a fictitious dependence 
on them). 



Theorem 3.1 The twisted u-Cauchy kernels of Sect. 3.1 are given by 



where 

Fr,,n{p, q) 



e 



{p 



G 



r„-(2^-l)/C + X)^^0^(^ 



(X - r„ - {2v - i)/c) ;A V 0a - cj ] 



-2iiTtj 



n 



9A(g-ooi+^) 

L k=l 0A(g- 00^ 



(3-69) 



(3-70) 

n > 

(3-71) 
n < 



77,... 



G C5 



^ = 2^tju(cj) - 2 



(3-72) 
(3-73) 



Correspondingly the dual {u, 1 — v)-differentials Pn,'^n 3re given by 



, , Tn{q)Y\U®i^{q-c,)-^'^'^ ^ 
Pn{q) = — 



c„eA((?-oo,^ ^ 



n+li 



T^nip) 



Trr\p)Y{UQA{p-c,) 

0a(P - OO^n+l) 



2iTTti 



-G 



-A 
-B 

A 
B 



q + Tn - (21^ - 1)IC - X - 00 



n+1 J e ^™ ' 



/iA(g) 



p - r„ (2z^ - i);c X - 00 



{+) \ -2iTTf r, 

n+1 ) e ^™ ' 



/ia(p) 
(3-74) 
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e 


'a 

B 


((2z^- i)/c + x-r„)e 


'a 

B 


({2,/-i)K + i:-r„+i) 











(3-75) 



They are defined up to rescaling by X, X ^ respectively (corresponding to different cfioices of tfie 
common basepoint po in the integrations above) and satisfy 



res pnl^m 
cxd(+) 



(3-76) 



Proof. The proof is a straightforward check that the proposed expression satisfies the defining 
properties in Eqs. (3-54). The expression T„(g) in eq. 3-71 has the following properties; 

• it has a tensor-weight of — ^ in the variable q, i.e. can be written in a local coordinate as 
f{z){dzy~2 and some multivaluedness around nontrivial cycles; 

• it has zeroes at r„ and poles at X; 

• it is projectively-independent (i.e. independent up to a multiplicative constant) of 
as long as they are not chosen so that the divisors appearing in the numerator are special. 
Indeed two different choices give functions with the same multivaluedness and the same divisor 
properties, hence proportional by a constant. 



r"(pj E(p,q) tensor weight u \x\ q and 1 — v \n p and the divisor 



The expression Gn{p,q) 
properties 

(Gn(p,*)) >r„-X-p , (G„(.,g)) > -r„ + X-(? (3-77) 

where r„ are defined in eq. (3-53). The remaining pieces of the formula make the final expression 
single valued and with the correct essential singularities. Q.E.D. 

As a corollary we derive the promised relation between the spectral bidifferential and the Szego 
kernels 

Corollary 3.3 The spectral bidifferential defined in Eq. 3-49 is given by 



Sip, 



Rn{P,q)^n{q,P) 



G 



{p-q-en)Q „ (q-p-en) 



e 



{-enfE^ip,q) 



X-Tn- i2u- 1)/C 



In [15], p. 26 we find that 



3 ^2 J 

S{p,q) = Q{p,q)+ 



j,k=i 



(3-78) 
(3-79) 

(3-80) 
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where Q is the normalized fundamental bidifferential (also known as Bergman kernel), such that 

s n = o. 

Jaj 

Moreover, S{p,q) is the square of the Szego kernel with complex characteristics; specifically, if 
p, € are the (half)-characteristics of e = F + {2v — 1)/C — X + A + tB 

e = 2p + 2re (3-81) 

then the above can also be rewritten as 

, e[^](p-g)e[^] (g-p) ^ I I , 

S{p,q) = e[^] {QfE^{p,q) = 5p,,(p,(?)5,,,(g,p) , (3-82) 

where the Szego kernel with characteristics is defined by 

4 Finite band recurrence relations and commuting (pseudo) dif- 
ference operators 

Let us partition the polar divisor of X into two disjoint subdivisors X = X(+) + X(-) of degree d 
and R + 1 — d respectively : 

(X)_ = - ^ oo5+) - ooj"^ =■ (4-1) 

i=0 3=0 

deg(X) = R+l . (4-2) 

We will choose the divisor 

^(+) 

as our dualization divisor. 
This means that according to the general scheme in Sect. 3.4 we will choose all the points oo^^'' 
within (h ence the same symbol is used). Correspondingly, all the points oo^- ^ are chosen within 

In this expression the points ooj^^ are not supposed to be necessarily distinct (within the same 
subset), so that we can consider poles of arbitrary order for X. 
We will postulate that 

ooj.;), = ooj"^ , r:=R-d+l (4-3) 
ooj+i = oo5+) (4-4) 

and assume that X has at least two distinct poles (r > 1); the modifications for the case of a single 
pole are left to the reader. We also fix a third kind differential rj as in Sect. 3.1.1. 
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Define the divisors 



iln := < 



n+1 



-r + X - J]] 00}+^ + ^ oor' n > 











-r + X + - Yl ^ < ° 

j=n+l £=n 

deg{iin) = -{2u-l){g-l)-l 



5J„ := <^ 



n+l 



00 



{-) 



n > 



£=1 




,{-) 



T-X- ooy> + ^ 00^ ^ n < 

j=—n £=—n+l 

deg(53„,) = (2z.-l)(g-l)-l 



(4-5) 
(4-6) 
(4-7) 
(4-8) 



In the formulas above it is understood that the sums are zero if the ranges are empty. The main 
point of these definitions is that 



(4-9) 



and these, in view of Sect. 3.4, are precisely the divisors characterizing (up to scalar multiplication)^ 

7r„ and p„ 

(vTn) > iln a (1 — z^)-differential 
(Pn) > '^^n a z^-differential 
These two sequences span vector spaces dual to each other under Serre's duality; 



(4-10) 



res TTnPm — ^mn 

X(+) 



Another crucial point that motivates the choice of twisting divisors is that now 

'C{pj}-d<j<R-d = 'Hu{To - 2X) 



and more generally 



'C{T^j}d~'R+n<j<n+d — 'Hl-u{—^n) 
'C{pj}n-d<j<n+R-d = 'Hui^n — 2X) 



(4-11) 



(4-12) 
(4-13) 

(4-14) 
(4-15) 



and hence we can conveniently choose them as components of the vectors (p^ used in the general 
construction. 



'^Thc essential singularities described in Sect. 3.1.1 are implied. 
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In this fashion, the vectors t/j^ and ip^ are windows of consecutive R+ \ elements within a pair 
of infinite vectors 



vr. 



n+d 



Pn-d 



Pn+R~d 



Definition 4.1 The vectors ^, ^ will 
be called the wave-vectors. The 
dual windows are defined by 

'■= [^n-R+d, ■■■ , VTn+d]* , 
fn '■= [Pn-d, ■ ■ ■ , Pn+R-d] , 

where the entries, depending on the 
base divisor F = Fq and the line bun- 
dle associated to the differential rj are 
given by the expressions in Thm. 3.1. 



Proposition 4.1 The wave vectors satisfy a finite 
band recurrence relation 



= $*X 



(4-16) 



X 



where the doubly infinite matrix X has a finite 
band structure of the form indicated on the right, 

and ajin) := res pn+jXTTn. 

Proof. The fact that ^',<1> solve transposed recurrence relations is immediate from the residue- 
pairing 




res '^(^ 



res = X 



(4-17) 



where 1 is the infinite identity matrix. The shape of the matrix X follows from inspection of the 
divisor properties of 7r„ and X7r„. Q.E.D. 

Keeping this in mind we can prove 

Proposition 4.2 The Christoffel-Darboux pairing is given by the non-zero (i? + 1) x (/? + 1) 
block in 

K„:=[n„,X], (4-18) 

where n„ = diag(. ..,...,1,0,...) is the projector up to n (i.e. the zero entries on the diagonal 
start at the entry (n -|- 1, n + 1) j. 

Proof. Using the definition of the matrix ]K„ 

(Xip) - Xiq))MP,Q) = <pi{q)KMp) (4-19) 
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we find that the entries of Kn are given by 

0^n)a,b = res res {X{p) - X{q))Rn(.P,q)Pa{p)n(.q) 

a = n — d, ...,n + R — d; b = n — R + d,...,n + d (4-20) 

Note that the order of the residues is irrelevant because the integrand is regular on the diagonal 
p = q. Using relations (3-67) we conclude that the nonzero entries are 



res XpaiTb b > n + 1 , a < n — 1 

X(+) 

- res XpaTTh b < n — 1 , a > n + 1 



(4-21) 



Explicitly, using the notation X„m = res ipn+kXipn = Oikin), we have (recall also that ak(n) = 
for k > d and k < d - R) 









ad{n - d) 


■ Oid{n) 


K = Kn = - 


-ad~R~i{n + I) ■■■ 


-a-i{n + 1) 
-ad^ji_i{n + d-R) 







;4-22) 



This concludes the proof. Q.E.D. 

Since all basis elements ip^,ip^ and the coefficients of K.^ have been expressed in the most direct 
way as Theta functions or residues thereof, we have achieved our goal of providing a completely 
explicit expression for the solution of the inverse spectral problem described by formula (3-23) and 
Section 3.2. 

Of course one should consider only one member of the sequence, without reference to the full 
sequence: so, for example, we can identify the matrix in (3-23) with the zeroth term. 

4.1 Lax and ladder matrices 

As it was explained in Sect. 3.4, we have a sequence of Lax matrices {An{x)]n&i, and an 
intertwining sequence of ladder matrices {Cn{x)}ne'L such that 



An+l{x) = Cn{x) ^An{x)Cn{x) 



(4-23) 
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The ladder matrices C„ are linear in x as follows from the explicit formula in Sect. 3.2.3. In this 
context the formula reads 

= res ^n.f)^nipY^n 

^^+,{p) = Cn{X{p))xljM ■ (4-25) 

The reader should check that they have companion-like form as in [5], where the coefficients of 
X appear in the nontrivial row (column) of the ladder matrix. 

The situation is very similar to the recurrence relation satisfied by orthogonal polynomials and 
generalization thereof [5]: in that case the recurrence relations are typically of Hessenberg form 
(d = 1 in our setting). 

However in the case of (generalized) orthogonal polynomials, the ladder matrices induce Schlesinger 
transformations for the associated Riemann-Hilbert problem, whereas here they have simply the 
meaning of an elementary twisting of the divisor. Note that the matrix representing multiplication by 
X in the infinite basis of the wave-functions (or, dually, <&) is a finite-band matrix by construction; 
on the other hand multiplication by Y does not result in a finite band matrix, namely the matrix 

Ynm ■= res YlTmPn (4-26) 

X(+) 

is in general a "full" doubly-infinite matrix; it has a finite number of nonzero supradiagonals (corre- 
sponding to the degree of the subdivisor of 2) supported at X^^^) but in general the part below the 
diagonal is not finite-band. 

Nonetheless the following matricial identity holds 

[X,Y] = (4-27) 

which makes sense entry-wise since, X being finite band, the commutator involves only a finite num- 
ber of terms. The commutativity clearly follows from the fact that the two matrices represent com- 
muting multiplication operators. Thus we are looking at a pair of commuting pseudo-difference 
operators, where X is a bona-fide difference operator, while Y is not. The only case in which Y is 
of finite band structure as well is if the divisor of poles of Y coincides with the polar divisor of X (al- 
though of different multiplicity in general). We remark however that in the case of pseudo-difference 
operators, even if the matrix Y has no obvious shape it has nonetheless a hidden rank condition. 

Indeed, under some additional genericity assumptions we can factor Y into the inverse of a 
lower-triangular matrix and a finite-band matrix. 

To see this, let us separate the poles of Y into the poles 2)^. that are also poles of X and the 
"other" poles 2)o of degrees dx and do respectively. 
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In a generic situation, we can find a linear combination 7i^„ G C{7r„, . . . , 7r„_^^} whose divisor 
exceeds 2)o and hence YTTn will be a linear combination of 7r„+j, for |j| < d^; the actual shape 
depends on how the poles are distributed into the divisors 

Setting 2)1?^^ = 2)a; fl X^^^ and letting m± be the respective degrees, we see that there is a 
lower-triangular matrix L with d„ subdiagonals and a finite-band matrix H with m+ supradiagonals 
and m_ subdiagonals such that 



(4-28) 



This implies that all the submatrices below the main diagonal have rank at most do and this is our 
"hidden" rank condition. 

Remark 4.1 A more general construction is needed in application to large degree asymptotics of certain biorthogonal 
polynomials (we will pursue this in a different publication [8]). In this case the twisting of the base divisor is performed 
in a different way so that both X and Y are pseudo difference operators subject to similar rank conditions. 



4.2 Riemann— Hilbert problems 

Although it is outside of the scope of this paper, we would like to explain what makes this investigation 
of potential relevance in a study of large degree asymptotics of (multi)-orthogonal polynomials. 

We start with the observation that the windows are eigenvectors of a matrix An{x) 

which depends only on the value x = X{p); since there are R + 1 = deg{X) other points 
Pi{x), . . . ,pji^i{x) € C (generically distinct) with the same X-projection, the evaluation at those 
points provides a basis of eigenvectors. The points (and so the eigenvectors) are distinct away from 
the ramification divisor of the map X : C ^ C and hence the sections pi : C C are well de- 
fined only on a suitable simply connected domain obtained by removing some cuts originating at the 
branch-points of the X-projection. 

The matrices 

V'n(Pl(2;)) „-2i7rf"<^'n i^niPR+li^)) 



Pl{x) :-- 



dXiplix)) 



dX(pR+i(x))l- 



-2i7r/ 



PI (a;) , 



dX(pi(x))- *^ 



(4-29) 



(4-30) 



AX(pn+^(x)Y 

are inverses of each other. They solve a Riemann-Hilbert problem with quasipermutation mon- 
odromies around the branch-points of X (due to the permutation of columns and to the multival- 
uedness of columns as functions on the spectral curve itself) and diagonal multivaluedness around 
the X-projection of the poles of rj (due to logarithmic singularities ). It would not be difficult, but 
too long, to spell out in detail the Riemann-Hilbert data, as they include some growth conditions at 
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infinity and at the branch-points. Note in general that we can expect singularities at a branch-point 
Xo of order k of the form (x — Xo)~ for Pr and of the form {x — Xo)~~ for Pl, or combination of 
singularities of this type if there are more than one ramification points on the spectral curve above 
the same branch-point. 

Riemann-Hilbert problems of this sort have been used in [10, 8] in the asymptotic analysis of 
certain (bi)orthogonal polynomials (for the case v = ^). 

The point of contact between the above RHP and the ones satisfied by (multi)orthogonal polyno- 
mials in the asymptotic regime is that such problems with quasi-permutation monodromies appear 
when the original RHP is "normalized" by the use of a suitable collection of (7-f unctions [12]. 

5 Commuting (pseudo)— difference operators in duality related 
to the two— matrix model 

We consider now a particular case which is of relevance for the asymptotic analysis of the biorthogonal 
polynomials for the so-called "two-matrix model" [18, 5]; we will remark later on what are the choices 
of the tensor weights and divisors which are more strictly relevant to that situation. 

The restriction will be that X and Y share the same polar divisor in the specific form [4] 

{X) > —oox — d2<X)y , iY) > —dioox — ooy . (5-1) 

We use cx)(+) = 00^;, cx)(~) = oOy and the same general framework used earlier, with a (generic) 
divisor F of degree {2iy — l){g — 1) + d2 and the third-kind differential rj. 

The wave-vectors are then characterized (up to constants) by the formulas in Thm. 3.1 suitably 
specialized; for reader's convenience we recall the divisor properties 

(vTn) > -T - noox -I- (n -I- d2)oOy , (p„) >T + {n- 1)qOx - {n + 1 + ^2)00^ (5-2) 

res TTnPm = ^mn ■ (5-3) 

Since the functions X,Y share the polar divisor, the sequence of wave-functions also satisfies a finite 
band Y-recurrence relation and hence the matrices X, Y can be thought of as two commuting 
difference operators[25]; denoting 

Ym,n ■■= -T^esYpmTTn (5-4) 

= -Y^ , y^>* = -^-^Y (5-5) 

we see that the two matrices X, Y are finite-band Hessenberg matrices with d2 subdiagonals for X 
and di supradiagonals for Y. 
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On the other hand we could have switched the roles of X and Y, used a (generic) divisor T of 
degree {2i/ — l){g — 1) + di and a differential rj and repeated the whole construction so as to get 
another pair of sequences of wave-functions 

(TTn) > -r - nooy + (n + di)QOx , (pn) > T + {n - l)oOy - (n + 1 + di)oOx (5-6) 

TeS7TnPm = Smn (5-7) 



Ynm = res YlTnPm , ^nm = " reS XTTnPra (5-8) 

OOy OOy 

It is clear that in general the matrices X, Y and their tilde-counterparts have nothing in common 
except the shape. In the formal asymptotics of biorthogonal polynomials the two matrices should 
however be the same [5]: we will see below that this implies certain constraints on the divisors r,r 
and the differentials r],rj. 
Suppose that 

r/-hry= dF = exact differential , T + f - X - Tl = {ly + i' - 1)C (5-9) 

where C is a canonical divisor. This means that there exists a {u + i' — l)-difFerential whose divisor is 
the one above. Also we assume that i], rj have the same polar divisor (in particular opposite residues). 
Under these conditions we have 

Proposition 5.1 If the divisorsT, F and differentials ri,rj are dual in the sense of eq. (5-9) then there 
exists a {u, v) -bidifferential £ which we call the Laplace kernel and a {I — v,l — T>)-bidifferential 
£ which we call the co-Laplace kernel with the properties 

(£(p, q))p >T -oox- d200y - q , {£{p, q))g >T- dioo^ - oOy - p (5-10) 
{£,{p, q))q > -T + oox + d200y - p , {£,{p, q))p > -T + dioox + oOy - q (5-11) 

Along the diagonal p = q they behave as 

£ ~ dz^d/^^^ + . . . (5-12) 
z — z' 

£ - dzi^'^dz''"^^^ + . . . (5-13) 
z — z' 

where lo := f{z)dz'^~^'^^^ and Q := f{z)dz^~'^^'^ are invariantly defined differentials of the indicated 
weights. 

Proof Similarly to the construction of the i/-Cauchy kernel, we split the divisors r,r into 

deg=q— 1 
da+l 2i/-l ^ ^ 

r=E7,+ E (5-14) 

j=l a=l 
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3=1 a=l 

and choose 2v — I points and 217—1 points ^j, (again, the formulae will depend only projectively 
on those). We rewrite formulae (3-70, 3-71, 3-74) specializing them twice: 

• the first time with oo^^ = oox, oon ^ = oOy,X = oox -|- d20Oy\ 

• the second time with oo^i^^ = ooy, ooi ^ = oOa;,2) = ooy + dioox- 
Duality (5-9) implies that 

A=-A, B = -B ,tj = -tj ,u{T + {2v-l)K-?i)=u[-f-{2V-l)K + 'X)^ (5-16) 
so that 

Fr,fl{p,q)=Frjfl{q,p) (5-17) 

with the quantities being defined as in eq. (3-70) and the above specializations. Then a direct 
inspection shows that (Tq defined as part of (3-71)) 

2i7rF(p)-2i7r r»7 



£(p,g) = ^^|^|^F^,o(p,<?)exp(-2^^ |% + 2.^F(g)j (5-18) 

£(p, q) = a^^lM^ e^- (5-19) 

T,{p)To{q)E{p,q) 

The {v + v — l)-difFerential and {\ — v — z?)-differential advocated for in the proposition are (up to 
multiplicative constant) 

uip) := ro(p)ro(p)e2-^(f) , Qip) = 4t , (5-20) 

where the fact that these expressions are single-valued follows once more from (5-9). Q.E.D. 

Corollary 5.1 The normalizations of the four sequences of wave-functions can be chosen in such a 
way that 

res £(C, q)TTn{0 = ^{q)'^n{q) = Pn{q) (5-21) 

res £(e,p)p„(0 = S(p)Pn(p)=^n(p) (5-22) 

res 2.{p,i)TTn{i) = ^{p)T^n{p) = Pn{p) (5-23) 

res 2.{q,i)pn{i)=^{p)pn{(l)=T^n{q) (5-24) 
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Proof. We first note that the differentials appearing in all the residues above are meromorphic 
since the essential singularities (by construction) cancel out. The only poles are a priori at cxo^. ^ and 
hence the sum over the two residues is (minus) the residue along the diagonal. Thus the statement 
of the theorem, keeping into account that ujuj = 1, amounts to checking that 

7r„(p) = Uj{p)pn{p) , Pn{p) = L0{p)TTn{p) ■ (5-25) 

The check is a straightforward computation using the explicit expressions (5-20) and the expressions 
for the two dual sequences derived from the suitable specializations of eqs. (3-74): the tilded 
sequences defined in (5-25) coincide with those obtained by specialization of (3-74) up to a rescaling 
Pn ^nPn, TTn ^ x~'^" (which leaves invariant the duality (3-76)). Q.E.D. 

Corollary 5.2 If the duality (5-9) is satisfied then the matrices representing multiplication by X and 
Y are the same (up to a transposition and a sign) in the two dual bases. 

Proof. Indeed 

X„m = res X-KnPm = reS Xuj-Kni^Pm = reS XpnTTm = - reS XpnTTm = -^mn (5-26) 

OOx OOx OOx OOy 

Y„,m = res YTTnPm = reS YaJTrnLOpm = res YpnTTm = - res YpnTTm = - Ymn (5-27) 

Q.E.D. 

There are two Lax-matrices (see [5] for the analysis for biorthogonal polynomials): An{x) of size 
{d2 -I- 1) X {d2 + 1) and Bn{y) of size (di -1- 1) x (di -I- 1) which share the characteristic polynomial 
(by construction). 

In addition, due to Corollary 3.3 and the duality 5-9 the two spectral bidifferential coincide. 

This means that if we denote by A{x) the Lax matrix reconstructed using X{p) as spectral 
parameter (of dimension d2 + 1) and B{y) the {di + 1)^ Lax matrix obtained by using instead Y{p) 
as spectral parameter and X{p) as eigenvalue, we have the identity 

det(xlrf^+i - B{y)) = cdet(ylrf,+i - A{x)) 

dxdx' Tr (^(A-?/)(x)(^-y)(x')) 
{x - x'Y (^(A:ry){^jc)^ TV ((l^)(xO) 
dydy' TT({B^x){y){^x){y')) 
(y - y'? T\- [{B^x){y)^ Tr [{B^x){y')^ 
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(5-28) 



(5-29) 



Conclusion. We conclude the section by pointing out the specialization that will be of use in the 
study of the asymptotics for biorthogonal polynomials appearing in the two-matrix model; indeed, 
the choice of the differentials 7],rj above was too generic. 

The relevant case would be u = u = in this case F, F have degrees d2,di respectively and one 
should choose them as F = ^200^, F = dioo^- The differential i] is then defined by linrj = NYdX 
and the dual one by 2i7rrj = NXdY so that in (5-9), liirF = NXY. 

The parameter N in the biorthogonal polynomial context is a large parameter (corresponding to 
the size of the underlying matrix model) and the degree of the polynomials whose asymptotics we 
are interested in, is n = + r with r € Z arbitrary but not scaling with N , i.e. bounded. 

In the spirit of this paper we should think of as a deformation of the line bundle rj while the 
"fluctuations around the Fermi level" (using a common analogy in the physically oriented literature) 
should be identified with r and give rise to the generalized Toda lattice (in this case it is the 2-Toda 
lattice). 
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